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I. INTRODUCTION 



III. EVENT SELECTION 



Experimental data for the processes tpi^S) 
Tr^J/tjj, vJ/i'j and 7Xci,2 are scarce and were mainly 
taken in the 1970s and 80s [1-6]. The branching frac- 
tions from different experiments do not agree well, 
and the ir'^J/ip channel is measured with low preci- 
sion. In particular, improved branching fractions for 
i/'(2S') — > 7Xcj are very important for the measure- 
ment of XcJ decay branching fractions using ipi'^S) 
data. All these appeal for high statistics measure- 
ments of these channels. 

In this paper, we report on the analysis of tp{2S) 
-K^J/ip, rjj/tjj and 7Xci,2 decays based on a sam- 
ple of 14.0 X 10*^ ip{'2S) events coUcctcd with the 
BESII detector. The first two decays are impor- 
tant to test various theoretical predictions for the ra- 
tios - r(^(2S)^7r°J/V.) _ r(T'^^T) , 

R" = -rt^Usr'^^i] I ^ ■ The ratio R has been calculated 
by different theoretical approaches [7], [8] [9], and the 
ratios R' and R" have been predicted in the framework 
of the QCD multipole expansion mechanism [10] [11]. 



II. THE BES DETECTOR 

The Beijing Spectrometer (BES) detector is a con- 
ventional solenoidal magnet detector that is described 
in detail in Ref. [12]; BESII is the upgraded ver- 
sion of the BES detector [13]. A 12-layer vertex 
chamber (VC) surrounding the beam pipe provides 
trigger information. A forty-layer main drift cham- 
ber (MDC), located radially outside the VC, provides 
trajectory and energy loss {dE/dx) information for 
charged tracks over 85% of the total solid angle. The 
momentum resolution is (7p/p = 0.017-\/l + {p in 
GeV/c), and the dE/dx resolution for hadron tracks 
is ~ 8%. An array of 48 scintillation counters sur- 
rounding the MDC measures the time-of-flight (TOF) 
of charged tracks with a resolution of ~ 200 ps for 
hadrons. Radially outside the TOF system is a 12 ra- 
diation length, lead-gas barrel shower counter (BSC). 
This measures the energies of electrons and photons 
over ^ 80% of the total solid angle with an energy 
resolution of ue/E = 22%/y/E {E in GeV). Outside 
of the solenoidal coil, which provides a 0.4 Tesla mag- 
netic field over the tracking volume, is an iron fiux 
return that is instrumented with three double layers 
of counters that identify muons of momentum greater 
than 0.5 GeV/c. 

A GEANT3 based Monte Carlo program with de- 
tailed consideration of the detector performance (such 
as dead electronic channels) is used to simulate the 
BESII detector. The consistency between data and 
Monte Carlo has been carefully checked in many high 
purity physics channels, and the agreement is quite 
reasonable. 



The data sample used for this analysis consists 
of (14.00 ± 0.56) X lO'^ V(2.S) events [14] collected 
with the BESII detector at the center-of-mass energy 
-y/s = M^(^2S)- The channels investigated are tp{2S) 
decaying into [tt'^ , rj) J / ip and 7X01,2, with Xci,2 de- 
caying to -yJ/ip, tt" and t] to two photons, and J/tp to 
lepton pairs. They all have a final 77/+/" {I = e,/i) 
event topology. 



A. General selection for 77/"*"/ events 

A neutral cluster is considered to be a photon can- 
didate if it is located within the BSC fiducial region 
(|cos0| < 0.75), the energy deposited in the BSC 
is greater than 50 MeV, the first hit appears in the 
first 6 radiation lengths, the angle between the clus- 
ter and the nearest charged track is greater than 14°. 
Each charged track is required to be well fit by a 
three-dimensional helix and to have a polar angle, 0, 
within the fiducial region | cos < 0.8. To ensure 
tracks originate from the interaction region, we re- 
quire V^y = ^Jvf+V^ < 2 cm and \V^\ < 20 cm, 

where V^, Vy, and Vz are the x, y, and z coordinates 
of the point of closest approach of each track to the 
beam axis. 

Events with two charged tracks and two or three 
photon candidates are subject to further selection cri- 
teria. The two charged tracks are identified as an 
electron pair or muon pair if < 5 < 0.6 or S > 0.9, 
respectively, where 



,E, 



scl 



Pi 



E. 



sc2 



P2 



1)2 



and p and Egc are the momentum and the deposited 
energy in the BSC of a charged track. 

A five constraint (5C) kinematic fit to the hypoth- 
esis iJj(2S) 77?"^Z~ with the invariant mass of the 
lepton pair constrained to the J/ip mass is performed, 
and the fit probability is required to be greater than 
0.01. For events with three photon candidates, the 
combination of two photons having the smaller x^ is 
chosen. Fig. 1 shows a scatter plot of the invariant 
mass of the reconstructed J/tp and the photon with 
higher energy {M^^^j/^) versus the invariant mass of 
two photons (Mj^) for the 77/x+/x~ final state. The 
corresponding plot for the 77e+e~ state is very sim- 
ilar. The 77, Xci, and Xc2 signals are qiute promi- 
nent, while the tt*^ signal is much less so. The cor- 
responding plot for 200,000 Monte Carlo simulated 
'>p{2S) — > TT^'jT^J/'tjj events, which is the main back- 
ground for the studied channels, is shown in Fig. 2. 
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(5.1% for 77e+e and 4.3% for 77/1+^ ) are included 




FIG. 1: M^i^^j/^ versus M-y^ after general selection of 
77/i^/i~ events. 



FIG. 3; Two photon invariant mass distribution for can- 
didate Tpi^S) Tv'^J/'ip events for (a) 77e~'"e~ and (b) 
77^t+p". 
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FIG. 2; M^i^^j/^ versus M^., after general selection for 
200,000 ^l){2S) TT°n°J/'4> Monte Carlo events (77/i+/i" 
final state). 



B. Selection of ip{2S) ■K^J/'ip 

To remove the huge background from ip{2S) 
7Xci,c2 under the il^{2S) 7r° J/-0 signal, we require 
M^^ j/^ to be less than 3.49 or greater than 3.58 
GeV/c^. Fig. 3 shows, after this requirement, the dis- 
tribution of invariant mass, Mj^, where the smooth 
background is contributed by ^{2S) 7Xci,2 and 
-0(25*) 7r°7r°J/-0. A Breit Wigner with a double 
Gaussian mass resolution function to describe the 7r° 
resonance plus a third-order background polynomial 
is fitted to the data. The fit gives N^°^- = 123 ± 18 
for the 77e+e~ state and iV7+, - = 155 =t 20 for 
the 77/i+/i^ state. In the fit, the mass resolution 
and the area ratio of the two Gaussians are fixed to 
the values determined by the Monte Carlo. The fit 
is also performed with a background function deter- 
mined by Monte Carlo simulated ip{2S) — > 
and V'(25') 7Xci,2 events that satisfy the selection 
criteria, where the two processes are weighted accord- 
ing to their branching fractions. The differences in the 
number of events obtained with the two backgrounds 



C. Selection of 1^(23) — >■ rjJ/i' 

In this channel the main backgrounds are from 
ipi^S) ~* TT^TT^J/ip and 7Xci.c2- By requiring 
M^^^j/,^ < 3.49 GeV/c^, most background from 
ip{2S) 7Xci,c2 is removed. The resultant plot, 
shown in Fig. 4, shows a clear 77 signal superimposed 
on background, mainly from '4'(2S) — s- tt'^tt'^ J/t/j. A 
fit is made using a Breit- Wigner resonance convo- 
luted with a mass resolution function for the rj sig- 
nal plus a polynomial background, where the width 
of the rj is fixed to its Particle Data Group (PDG) 
value [15] and the background function is deter- 
mined from ip{2S) — > 7r°7r° J/i/; Monte Carlo simulated 
events that satisfy the same criteria as the data. The 
fit yields N^+^- = 2465 ± 101 for the 776+6" state 
and iV^+^- ='^3290 ± 148 for the 77/1+ A*" state. The 
fitted values of the 1] mass are 547.6 ± 0.1 MeV/c^ 
for the 776+e~ channel and 547.7 ± 0.1 MeV/c^ for 
the 77/i+/i" channel, consistent with the PDG value 
within 2(7. 

A fit using a fourth-order background polynomial 
with parameters free is also performed to estimate the 
systematic error due to the background shape. This 
error is negligibly small. 

D. Selection of ■ip{2S) — > 7Xci,c2 

The processes tp{2S) 7r° J/i/j, ijj/^p, and n^TT^J/ip 
contribute to the background for this channel. By 
requiring M^^ < 0.53 GeV/c^, most of the back- 
ground from 'ip(2S) rjJ/ip and a significant por- 
tion from ip{2S) tt^tt^J/iP are rejected. Fig. 5 
shows the M^^^j/^ distribution for candidate ip{2S) 
7Xci,c2 events. The remaining background is mainly 
due to ip{2S) — > TT^ Ti^ J / ip . The contribution from 
ip{2S) — > 'K^J/i) is negligible due to its tiny branch- 
ing fraction. Fig. 6 shows the M^^jj^ distribution for 
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FIG. 4: Two photon invariant mass distribution for can- 
didate il}{2S) vJ/i^ events for (a) 776^6^ and (b) 
77/i+^i". 

iIj{2S) ■K^n^J/ij} Monte Carlo simulated events be- 
fore and after the M^^ < 0.53 GeV/c^ requirement, 
the latter one is taken as the background shape in 
the fit. Two Breit-Wigner resonances convoluted with 
mass resolution functions plus a background function 
are fitted to the data. The widths for the Xci,2 are 
fixed to the PDG values, and the mass resolution func- 
tions are determined by Monte Carlo simulation. The 
fit yields: 

iV\^i = 5263 ± 124, N^f = 2512 ± 82, 

e+e~ ' e+e ' 



where Y stands for the final state, X the intermediate 
state, and e the detection efficiency. The branching 
fraction of X ^ F is taken from the PDG. 



A. Detection efficiency 

The detection efficiency is the product of the trig- 
ger efficiency Ctrg and the reconstruction-selection ef- 
ficiency €rs- For the BES detector, the trigger ef- 
ficiency for hadronic events is 1.000 ± 0.005 [16]. 
The reconstruction-selection efficiency is determined 
by Monte Carlo simulation. For the signal chan- 
nels studied, generators taking into account phase 
space, angular distributions, and final state radiation 
are used for the event simulations. For the channel 
V'(2S') ■n^TT^J/tp, the common background for all 
signals, we use a generator, which gives the correct 
dipion mass and angular distributions [17]. 

For each of the channels analyzed, 50,000 Monte 
Carlo events are subjected to the same reconstruction 
and event selection as used for the data to determine 
the detection efficiencies, which are listed in Table III. 



B. Efficiency corrections and systematic errors 



- 6752 ±178, iV„Y^_ 



3358 ±96, 



with the fitted masses of Xci and Xc2 equal to 3510. 9± 
1.0 MeV/c2 and 3555.9 ± 1.0 MeV/c^, respectively, 
consistent with the PDG values. 




3.6 

.(GeV) 



Because the Monte Carlo is imperfect, it is neces- 
sary to correct the detection efficiencies obtained from 
simulations for the differences between the Monte 
Carlo simulation and the data. Differences come from 
the efficiencies of MDC tracking, particle identifica- 
tion, photon identification, and kinematic fitting. In 
addition, the uncertainties of the background shapes 
(estimated in Section 3), the number of tp{2S) events, 
and the branching fractions of the intermediate states 
also contribute to the final systematic error. 

To investigate the difference in the lepton track 
efficiencies of the Monte Carlo simulation and the 
data, the lepton pair sample from the decay ^(2S) 
n^TT" Jf-ip, J/ip , which closely simulates the 

behavior of the lepton pair in the channels under 
study, is used. This study finds the tracking efficiency 



FIG. 5: Invariant mass j/^ distribution for candidate 
iIj{2S) 7Xci,c2 events for (a) 77e'''e~ and (b) 77/^^/1". 



IV. BRANCHING FRACTION 
DETERMINATION 



For iP{2S) 
mined from 



X, the branching fraction is deter- 




B(V'(2S') ^ X) 



•0(25) 



B{X Y) ■ e{ip{2S) ->X~>Yy 



FIG. 6: Invariant mass M. 



distribution for Monte 



Carlo simulated ip{2S) — > J/tp events (77/i^^ final 
state), (a) Before the Mjj < 0.53 GeV/c^ requirement, 
(b) After the M-y-y < 0.53 GeV/c^ requirement. 
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correction factor is 1.012 ± 0.009 for e+e" pairs and 
1.002 ± 0.008 for ^"""/x" pairs. For charged particle 
identification, S is used to separate e+e" and lJ.'^li~ 
pairs. The same lepton pair sample is used to de- 
termine the particle identification efficiency difference 
between Monte Carlo simulation and data by deter- 
mining efRciences for each with and without this par- 
ticle identification requirement. The correction fac- 
tor is found to be 0.951 ± 0.008 for e+e~ pairs and 
0.972 ± 0.006 for fi+ fjr pairs. 

For the photon selection used, studies show that the 
efficiency difference between data and Monte Carlo is 
2% for each photon [18]. We take this difference as 
the systematic error in photon selection, and no cor- 
rection to the efficiency is made. In addition, the rib 
in the BSC causes an inefficiency in photon detection. 
The systematic error due to the rib efficiency, listed 
in table I, is obtained by comparing results with pho- 
tons in the rib region removed with those when they 
are not removed. 

The systematic error due to kinematic fitting comes 
from the differences between data and Monte Carlo 
simulation in the measurements of track momen- 
tum, the track fitting error matrix, and the pho- 
ton energy and direction. For the charged track 
part, the difference is estimated using the '(p{2S) 
tt'^tt" J/xp, J/ip sample. For the photon part, 

a careful calibration of the neutral cluster information 
in the BSC is performed and the difference with and 
without the calibration applied to both the data and 
Monte Carlo is used to determine the systematic error 
in this part [19]. 

Table I summarizes the efficiency correction fac- 
tors and uncertainties from all sources, while Ta- 
ble II lists the systematic errors for the channels under 
study. The branching fractions and corresponding er- 
rors for all intermediate state decays are taken from 
the PDG [20]. 



the Crystal Ball results [6]. Much of the systematic 
error on B{xp{2S) — > 7Xci,c2) comes from the uncer- 
tainties on B{xci — * and B{xc2 jJ/i')- 

Using Partially Conserved Axial-vector Currents 
(PCAC), Miller et al [8] predicts: 



R 



(1) 



wherer = (ma— m„)/(ms — 0.5-(md+m„)) andp,r and 
Pr, are the tt and rj momenta in the il>{2S) rest frame. 
With the conventionally accepted values of TOs = 150 
MeV/c2, nid = 7.5 MeV/c^, m„ = 4.2 MeV/c^ given 
by Weinberg [21], the ratio R equals 0.0162, which is 
smaller than our measurement (0.048 ± 0.005) [22] by 
a factor of three. Based on an effective Lagrangian 
approach, Casalbuoni et al [9] obtain an improved ex- 
pression 



R='^(i^ 



16 \pn 



1 



2B 



A/, 



3A 



1 -I- ^ 



B Xf, 

— : 



(2) 



in which A characterizes the rj—rj' mixing, B/A is a not 
yet determined parameter in the effective Lagrangian. 
,U = (130 ± 5) MeV is obtained from PDG. Using the 
approximation [23] 




rris - 



tan^ 



(3) 



where Op « 
we obtain A ; 



-20° [15] is the rj — rj' mixing angle, 
1.91 GeV. With our measured value of 



R, we infer the parameter ^ equals — 1.42 ± 0.12 or 
-3.11 ± 0.15 in Equation ( 2). 

In terms of QCD multipole expansion, Kuang 
et al [11] predict the ratio 



C. Results and discussion 

Using the fitting results and the efficiencies and 

correction factors for each channel, we determine the 
branching fractions listed in Table III. We also obtain 
the product branching fractions 

B(V(25) ^ 7Xci) • Bixci - iJ/^) = 
(2.81 ±0.05 ±0.23)%, 

B{^{2S) ^ 7Xc2) • B{xc2 ^ lJ/4>) = 
(1.62 ±0.04 ±0.12)% . 

Our B{tp{2S) n^J/tp) measurement has im- 
proved precision by more than a factor of two com- 
pared with other experiments, and the BES ip{2S) 
rjJ/il) branching fraction is the most accurate single 
measurement. Our B{il;{2S) TT°J/ip) agrees better 
with the Mark-II result [5] than with the Crystal Ball 
result [6], while B{xp{2S) 7Xci,c2) agrees well with 



~ \mj • \p,ms))) ■ \fms))) ' ^ ^ 

where f{ip{2S)),f{T'), and /(T") are the transition 
amplitudes of ^{2S) J/ijjin:^ T' — s- Ttttt, and 
T" Ttttt, respectively, which depend on the po- 
tential model describing the heavy quarkonia. Taking 
the QCD motivated Buchmitller-Grunberg-Tye poten- 
tial [24] as an example, the predicted values are 
R'^^rj, = 0.0025 and i^'^^j, = 0.0013. With our mea- 
surements of B{tp{2S) V^/i^) and PDG values of 
r(V'(2S')), r(T' ^ Try) and r(T" ^ T77), we obtain 
R'^^p < 0.0098 and i?"^^ < 0.0065, which are consis- 
tent with the predictions of Equations ( 4) and ( 5). 

The BES collaboration acknowledges the staff of 
BEPC for their hard efforts. The authors also thank 
Prof. Y. P. Kuang for enlightening discussions. This 



6 



TABLE I: Efficiency correction factors. 



Channel 




rtJ/ip 


Final state 


77e~''e 7TA*~'~M 


77e~'~e y^fx^fi 


Track selection 
Particle ID 

5-C fit 
7 efficiency 
r>oU riD 


1.012 ±0.009 1.002 ±0.008 

0. 951 ± 0.008 0.972 ±0.006 
1.000 ±0.014 1.00 ±0.02 

1.00 ± 0.04 1.00 ±0.04 

1. UUU IE U.Uzo l.UUU ±: U.Uo4 


1.012 ±0.009 1.002 ±0.008 

0. 951 ± 0.008 0.972 ± 0.006 
1.000 ±0.016 1.000 ±0.038 

1.00 ±0.04 1.00 ±0.04 

1. UUU ± U.Uoi i.UUU it U.UoD 


Total correction 


u.yoz ^ u.uou u.y/4 ^ u.ud/ 


u.yoz ^ u.uoo u.y/^i: ^ u.ud* 


Channel 


7Xci 


7Xc2 


Final state 


77e"'"e~ 77/x"'"/j~ 


77e"'"e~ 77/^"*"/^" 


Track selection 
Particle ID 
5-C fft 

7 efficiency 
BSC rib 


1.012 ± 0.009 1.002 ± 0.008 
0.951 ± 0.008 0.972 ± 0.006 
1.000 ± 0.015 1.000 ± 0.049 

1.00 ± 0.04 1.00 ±0.04 
1.000 ± 0.043 1.000 ± 0.040 


1.012 ± 0.009 1.002 ± 0.008 
0.951 ± 0.008 0.972 ± 0.006 
1.000 ±0.018 1.00 ±0.052 
1.00 ±0.04 1.00 ±0.04 
1.000 ± 0.019 1.000 ± 0.024 


Total correction 


0.962 ± 0.061 0.974 ± 0.075 


0.962 ± 0.049 0.974 ± 0.070 



TABLE II: Systematic errors (%) 



Channel 




rjj/d) 


Final state 


77e"'"e 


77/i+/i 


77e+e 


77/i+/j 


efficiency correction 


6.3 


5.9 


5.7 


6.9 


Number of ip{2S) events 


4 


4 


4 


4 


B(7r°,»7 ^ 77) 


negligible negligible 


0.65 


0.65 


BiJ/ip e+e-,n+n-) 


1.7 


1.7 


1.7 


1.7 


background shape 


5.1 


4.3 


negligible negligible 


Total systematic error {%) 


l).2() 


cS.:)0 


7.20 


8.18 


Channel 


7Xci 


7Xc2 


Final state 


77e"''e 


77^i+^ 


776^6 


77/i+/i 


efficiency correction 


6.3 


7.7 


5.1 


7.2 


Number of 'ip{2S) events 


4 


4 


4 


4 


BiXcJ ^ iJ/i') 


8.5 


8.5 


8.9 


8.9 


B{J/i! e+e-,fi+fi-) 


1.7 


1.7 


1.7 


1.7 


Total systematic error (%) 


11.44 


12.26 


11.14 


12.25 



TABLE III: Results. Note that much of the systematic error on B{tj){2S) — » 7Xci,c2) is due to the uncertainty on 
B(Xcl,c2 7-^/V')- 



Cliaiiiicl 


7r"J/V' 




riJ/il) 


Final state 


77e+e~ 77/U+jU~ 


776+ e 




Number of events 
EfRciency(%) 
Sys. error (%) 

Correction factor 


123 ± 18 155 ± 20 
11.21 13.34 
9.68 8.77 
0.962 0.974 


2465 ± 101 
26.94 
8.54 
0.962 


3290 ± 148 
34.07 
8.40 
0.974 


BR (%) 


0.139 ± 0.020 ± 0.013 0.147 ± 0.019 ± 0.013 


2.91 ± 0.12 ± 0.21 3.06 ± 0.14 ± 0.25 


Combine BR (%) 
PDG (%) 


0.143 ± 0.014 ± 0.013 
0.096 ± 0.021 


2.98 ±0.09 ±0.23 
3.13 ±0.21 


Cliaiinel 


7Xcl 




7Xc2 


Final state 


77e"'"e~ 77^^^^ 


776+ e 


77^t+^- 


Number of events 
Efficiency(%) 
Sys. error (%) 

Correction factor 


5263 ± 124 6752 ± 178 
23.88 29.24 
12.23 12.45 
0.962 0.974 


2512 ± 82 
19.70 
12.10 
0.962 


3358 ± 96 
25.54 
12.44 
0.974 


BR (%) 
Combine BR (%) 
PDG (%) 


8.73 ± 0.21 ± 1.00 9.11 ±0.24 ± 1.12 
8.90 ± 0.16 ± 1.05 
8.4 ±0.6 


7.90 ± 0.26 ± 0.88 8.12 ± 0.23 ± 0.99 
8.02 ±0.17 ±0.94 
6.4 ±0.6 
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